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Abstract 


Investigation of annular obstructed azimuthally polarized Bessel-modulated Gaussian beam 
(QBG) in the focal region of high NA lens, based on vectorial diffraction theory. The numerical results 
show that the intensity distribution in focal region of the incident beam can be altered considerably by 
changing beam parameter (u) and introducing annular apodization (6). Beam parameter induces the focal 
splitting in transverse direction, while annular apodization leads to change in focal pattern along optical 
axis of the focusing system. More interesting, the focal splitting may be in continuous in certain case of 
incident beam propagating through aligned optical system which is suitable for application such as optical 


manipulation and optical trapping. 
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1. INTRODUCTION 


The laser beams with cylindrical symmetry in 
polarization have attracted very much for their 
interesting properties and applications called 
cylindrical vector beams. This kinds of beams can be 
generated by using active and passive methods (Zhan 
et al. 2009; Youngworth et al. 2000; Gao et al. 2007; 
Zhou et al. 2007). Due to the polarization symmetry, 
the electric field at the focus has unique polarization 
properties. K.S. Youngworth and T.G. Brown 
calculated cylindrical-vector fields (Youngworth et al. 
2000), which shows that, in the particular case of a 
tightly focused radially polarized beam, in the focal 
region the polarization shows large inhomogeneities, 
while the azimuthally polarized beam has purely 
transverse field in focal region. Focus shaping 
technique is also reported by using generalized 
cylindrical vector beams (Gao et al. 2007), in which a 
generalized cylindrical vector beam can _ be 
decomposed into radially polarized and azimuthally 
polarized components. And a generalized cylindrical 
beam can be generated from a radially polarized or an 
azimuthally polarized light using a two-half-wave-plate 
polarization rotator. Recently, the tight focusing 
properties of a polarized beam and a circularly, radially, 
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azimuthally or linearly polarized vortex beam by a high 
numerical-aperture (NA) lens have been discussed 
(Quabis et al. 2001; Helseth et al. 2004; Grosjean et al. 
2007; Ganic et al. 2003; Zhan et al. 2002; Helseth et al. 
2006; Zhan et al. 2006). The study of coherent 
superposition of radially and azimuthally polarized 
components has also generated interest (Bokor et al. 
2007). The polarization propagation and light intensity 
distribution in focal region plays an important role 
in many optical systems. It was demonstrated both 
theoretically and experimentally (Jia et al. 2005, Suresh 
et al. 2013). Optical intensity distributions in focal 
region have found applications in optical micro 
manipulation domain (Ashkin et al. 2003; Paterson et 
al. 2001; Gao et al. 2004; Caron et al. 1999; Hricha et 
al. 2005). Pure phase-shifting apodizer effect in focal 
region is reported (Gao et al. 2009). Nowadays, the 
Bessel-modulated Gaussian beams with quadratic 
radial dependence (QBG beam) have novel focusing 
properties in cylindrical coordinate system. Caron and 
Potvliege introduced QBG (Caron et al. 1999; Hricha 
et al. 2005; Wang et al. 2002; Belafhal et al. 2000; Lii 
et al. 2002), it has attracted much attention (Wang et al. 
2003). The axisymmetric QBG beam are in zeroth- 
order, which is usually referred to as the axisymmetric 
QBG beam (Wang et al. 2008), Young worth and 
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Brown calculated cylindrical-vector fields, near the 
focal region of an aplanatic lens, and briefly discussed 
some applications (Youngworth ef al. 2000). Zhan and 
Leger reported a focus shaping technique using 
generalized cylindrical vector beam. In our knowledge, 
almost all QBG beams in previous papers are in scalar 
form without considering the polarization distribution 
of optical field. In this paper we numerically study the 
focusing properties of azimuthally — polarized 
axisymmetric QBG beam with radial variance phase 
wavefront. Though the proposed method is simple and 
free from the complex interferometric methods used in 
the conventional techniques to shift the generated focal 
hole. The principle of the focusing radially polarized 
axisymmetric QBG beam with radial variance phase 
wavefront is given in Section 2. Section 3 shows the 
simulation results and discussions. The conclusions are 
summarized in Section 4. 


2. THEORY 


In the focusing system we investigated, 
focusing beam is azimuthally polarized axisymmetric 
QBG in radial (Youngworth et al. 2000; Zhan et al. 
2002; Dorn et al. 2003; Suresh et al. 2014) Therefore 
the cylindrical coordinate system (r,¢,0) the field 
distribution E(r,¢,0) of the azimuthally polarized 
axisymmetric QBG beam at the plane z = 0 is written 
as, 


E,(r,9,0) = E,(7,9,2= 0)-| cos(¢(r))-n, | (1) 


Where ¢(r) is the polarization angle from 
azimuthal direction, and is the function of azimuthal 
coordinate for the azimuthally polarized axisymmetric 
QBG beam, were n, is azimuthal unit vector. 


r r 
E,(r,.9,z=0)=Jy e ex| S (2) 
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Where Jo denotes the Bessel function of order 
Zero, Wo is the waist width of the Gaussian beam, i is a 
beam parameter which is complex valued in general. 
After some simple derivation (Tang et al. 2009), Eq. 2 
can be rewritten as 
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Parameter w = ,/f is called relative waist 
width, where r, is radius of incident optical aperture. 
NA is numerical aperture of the focusing system. @ is 
polar angle corresponding to radial coordinate. ¢(r) is 
the polarization characteristics of the focusing 
azimuthally polarized axisymmetric QBG beam and 
can be expressed as 


poco 2", ‘és 
tan (a) 

Where a = arcsin (NA/n) is the convergence 
angle of the focusing optical system. If the focusing 
optical system investigated is in air, i.e. refractive index 
n = 1. C is radial variance parameter indicating 
polarization radial degree the wave front phase 
distribution. Using the same analysis method as that in 
references (Youngworth et al. 2000; Gu et al. 2000) the 
electric field in focal region of azimuthally polarized 
axisymmetric QBG beam is, 


E(r, Q.z) — Ee, + Ee, Ee, 
(5) 


E,(7,9,2)= 2A] ,fcos(4) sin(9)E(A)¥, (krsin 0 )exp(ikz cos0)d0 
0 


E,(r,9,z)=E,(r,9,z) =0 
(6) 
Where A is a constant, J,(x) and J;(x) denote zero 
order and first order Bessel functions of the kind, k is 
the wave number & = 27/1, with i being the wavelength 
of illuminating beam. Were e;, ey and e, are unity 
vector. 

Substitute Eq. (3) and (4) into Eqs. (6), then 
substitute Eqs. (6) into Eq. (5), we can obtain the 
analytical total electric field in focal region. The 
optical intensity in focal region is proportional to the 
modulus square of Eq. (5), so the focusing properties of 
azimuthally polarized axisymmetric QBG beam can be 
investigated theoretically. 


3. RESULTS AND DISCUSSIONS 


Without loss of validity and generality, it is 
proposed that the numerical aperture of the focusing 
system NA = 0.95 and amplitude constant A = /. The 
focusing properties of azimuthally polarized 
axisymmetric QBG beam is investigated theoretically, 
for parameter C = 0 at different beam parameter (wu). 
Firstly, the intensity distributions in focal region of the 
azimuthally polarized QBG beam are calculated under 
condition of C = 0 and different uw, and are illustrated 
in Fig. 1. In this paper, it should be noted that the unit 
of coordinates in all figures in this article is k/. Focus 
refers to the maximum optical intensity peak, and focal 
shift denotes the movement of this focus. 


Fig. 1 shows the total intensity distribution of 
incident azimuthally polarized QBG in the focal region 
of high NA lens (NA = 0.95). It can be seen from Fig. 
1(a) in the focal region for uw = O, ie. the system 
generates two focal spot with equal intensities in radial 
direction forms focal hole of focal depth (1.4A) each. 
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The focus (focal hole) in the focal region splits along 
optical axis on increasing parameter uw, namely focal 
split phenomenon occurs and is shown in Fig. 1(b) for 
ue = 4. On increasing = 6 continuously, the focal 
pattern changes considerably and illustrated in Fig. 
1(c). There is two bright spot with two weaker spot on 
either side of bright spot in the focal region. We further 
increased u = 9, in order to study the effect of beam 


parameter pt and its evolution in the focal region is 
shown in fig. 1 (d). We observed that on increasing 
the generated focal hole splits along optical axis and 
forms four bright spot in the focal region, the bright 
spot becomes weaker and the weaker spot becomes 
stronger. 
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Fig. 1: Total intensity distribution of azimuthally polarized high NA lens (NA = 0.95) for (a) u=0, (b) u=4, 


(c) ¢=6 and (d) u=8 


Secondly, we analyzed numerically the effect 
of annular apodization for different parameter yw by 
keeping other parameter same as above. The effect of 
annular apodization on the incident radially polarized 
axisymmetric QBG beam is analyzed replacing E(@) by 
E(0)T(@) in Eq. (5). Where 7(0) is given by 


O fe 
rO)=1, a 


In order to study the effect of annular 
apodization we choose annular apodizer 6 = 0.85, by 
keeping other parameter same as above for different 
pupil beam parameter and it is shown in Fig 2. 
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Fig. 2: Total intensity distribution of azimuthally polarized high NA lens (NA = 0.95) with annular 
obstruction (6 = 0.85) (a) u=0, (b) u=3, (c) u=4, (d) u=6, (e) u=7, (f) u=8, (g) H=I and (h) w=10 
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Fig. 2 (a) shows the total intensity distribution 
contour plot of incident annular obstructed azimuthally 
polarized axisymmetric QBG beam when yu = 0, the 
focal region has only one focal hole (two focal spot in 
radial direction with FWHM (0.28) and its focal depth 
of FWHM (4.064) of each). Further increase in beam 
parameter u, results in focal splitting along axial 
direction at the focus and it is shown in Fig. 2 (b) for uw 
= 3, for u = 6, generates focal hole in the focal region, 
it is shown in Fig 2 (c). Fig 2 (d) for uw = 6, Fig. 2 (d) 
for u = 7, Fig. 2 (d) for u = 8, Fig. 2 (d) for uw = 9 and 
Fig. 2 (d) for u = 10. It is observed from Fig. 2 on 
increasing continuously, the focal pattern changes 
from one focal hole to two focal hole in the focal 
region. 


4. CONCLUSION 


Based on vector diffraction theory, tight 
focusing of a azimuthally polarized QBG beam by a 
high NA lens is studied theoretically. The total optical 
intensity distribution in the focal region is illustrated by 
numerical calculations without and with annular 
apodization for different beam parameter yu. The 
simulation result shows, on increasing beam parameter 
y the proposed system generates different focal pattern 
changes at the focus. We expect such a beam can be 
widely used in application such as optical trapping, data 
storage, biomedical imaging, laser drilling, and laser 
machining. 
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